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Abstract 

Energy detection is a widely used method of spectrum sensing in cognitive radio and Radio 
Detection And Ranging (RADAR) systems. This paper is devoted to the analytical evaluation of the 
performance of an energy detector over Weibull fading channels. This is a flexible fading model that has 
been shown capable of providing accurate characterization of multipath fading in, e.g., typical cellular 
radio frequency range of 800/900 MHz. A novel analytic expression for the corresponding average 
probability of detection is derived in a simple algebraic representation which renders it convenient 
to handle both analytically and numerically. As expected, the performance of the detector is highly 
dependent upon the severity of fading as even small variation of the fading parameters affect significantly 
the value of the average probability of detection. This appears to be particularly the case in severe 
fading conditions. The offered results are useful in evaluating the effect of multipath fading in energy 
detection-based cognitive radio communication systems and therefore they can be used in quantifying 
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the associated trade-offs between sensing performance and energy efficiency in cognitive radio networks. 


I. Introduction 

Detecting unknown signals is a critical topic in wireless communications and has attracted 
for years the attention of both Academia and Industry. The detection is typically realized in the 
form of spectrum sensing with energy detection (ED) constituting the most simple and popular 
method. The operating principle of ED is based on the deployment of a radiometer, which is 
a non-coherent detection device that measures the energy level of a received signal waveform 
over an observation time window. The obtained measure is subsequently compared with a pre¬ 
defined energy threshold which determines accordingly whether an unknown signal is present 
or not [[Il-l|3l. 

Detection of unknown signals over a flat band-limited Gaussian noise channel was firstly 
addressed by H. Urkowitz in dH. There, he comprehensively derived analytic expressions for the 
corresponding probability of detection, P^, and probability of false alarm, Pfa, performance 
metrics. These measures are generally based on the statistical assumption that the decision 
statistics follow the central chi-square and the non-central chi-square distribution, respectively. A 
few decades later, this problem was revisited by Kostylev [|5l who considered quasi-deterministic 
signals operating in fading conditions. 

Thanks to its low implementation complexity and no requirements for prior knowledge of 
the signal, the ED method has been widely associated with applications in RADAR systems 
while it has been also shown to be particularly useful in emerging wireless technologies such as 
ultra-wideband communications and cognitive radio flU, [|71. In the former, the energy detector is 
exploited for borrowing an idle channel from authorized users, whereas in the latter it identifies 
the presence or absence of a deterministic signal and decides whether a primary licensed user 
is in active or idle state, respectively. According to the corresponding decision, the secondary 
unlicensed user either remains silent, or proceeds in utilizing the unoccupied band until the 
primary user becomes again active [[D, |[8l. This opportunistic method has been shown to increase 
substantially the utilization of the already allocated radio spectrum, offering a considerable 
mitigation of the currently extensive spectrum scarcity ||9l- [IT5]l . 
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It is recalled that fading effects create a detrimental impact on the performance of wireless 
communication systems [fT^ - [|25ll and the references therein. Several studies have been devoted 
to the analysis of the performance of energy detection-based spectrum sensing for different 
communication and fading scenarios. Specifically, the authors in [|^ derived closed-form ex¬ 
pressions for the average probability of detection over Rayleigh, Rice and Nakagami-m fading 
channels for both single-channel and multi-channel scenarios. Similarly, the ED performance 
in the case of equal gain combining over Nakagami—m multipath fading has been reported by 
[l27]l whereas the performance in collaborative spectrum sensing and in relay-based cognitive 
radio networks has been investigated by [|^ - [[^ . A semi-analytic method for analyzing the 
performance of energy detection of unknown deterministic signals was reported in If^ and is 
based on the moment-generating function (MGF) method. This method was utilized in the case 
of maximal-ratio combining (MRC) in the presence of Rayleigh, Rice and Nakagami-m fading 
in [1^ as well as for the useful case of correlated Rayleigh and Rician fading channels in ll34]l . 
Finally, the detection of unknown signals in low signal-to-noise-ratio (SNR) over iT-distributed 
{K), generalized K {Kg) and the flexible rf—jj, and K—jj, generalized fading channels have been 
analyzed in [|35l-[|39l. 

Weibull statistical distribution was proposed in 1937 and became widely known by Walodi 
Weibull in 1951. It has been used in several fields of natural sciences and engineering in 
applications relating to reliability, failure data analysis, weather forecasting and modelling of 
signal level dispersion in RADAR systems. In the context of wireless communications, it has 
been employed in modelling multipath waves propagating in non-homogeneous environments. 
Weibull distribution is very flexible and has been shown to provide adequate fitting to empirical 
results from wireless channel measurements in both indoor and outdoor communication scenarios. 
Furthermore, it is capable of accounting for propagation in urban environments in cases, where 
Rayleigh distribution is inadequate, as well as in modelling path-loss in, e.g.. Digital Enhanced 
Cordless Telecommunications (DECT) system at 1.89 GHz Il40l - [|45]| . Il48l . Its flexibility is also 
evident from the fact that it includes as special cases the well-known Rayleigh and exponential 
distributions for a = 2 and a = 1, respectively. 

However, in spite of the undoubted usefulness of Weibull fading model, no studies related to 
the detection of unknown signals over such fading conditions have been reported in the open 
technical literature. Motivated by this, the present work is devoted to the analytical performance 
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evaluation of energy deteetion over narrowband Weibull fading ehannels. Speeifically, a novel 
analytie expression is firstly derived for the average probability of deteetion of the energy deteetor. 
This expression is rather simple both algebraieally and numerieally and is subsequently used 
to analyze the performanee of different fading conditions. As expected, the performance of 
the energy detector is highly dependent upon the value and variations of the involved fading 
parameters. This allows the quantification of the effect of fading in the system performance 
in various communication scenarios which can enable the determination of the required power 
levels for ensuring robust and accurate performance of energy detectors along with sufficient 
benefits, when possible, with respect to energy efficiency. 

The remainder of this paper is organized as follows: The system and channel model are 
described in Section II. The average detection probability of the energy detector over Weibull 
fading channels is analyzed in Section III and numerical results for each communication scenarios 
and discussions are provided in Section IV. Finally, closing remarks are given in Section V. 


II. System and Channel Model 


A. Energy Detection 

The received signal waveform in narrowband energy detection follows a binary hypothesis 
that can be expressed as |l34l eq. (1)], 



( 1 ) 




where s{t) is an unknown deterministic signal whereas h denotes the complex gain of the 


channel coefficient and n{t) is an additive white Gaussian noise (AWGN) process. The samples 


of n{t) are assumed to be zero-mean Gaussian random variables with variance NqW with W and 
Nq denoting the single-sided signal bandwidth and a single-sided noise power spectral density, 
respectively [l34ll . The hypotheses Hq and Hi refer to the cases that a signal is absent or present, 
respectively. The received signal is subject to filtering, squaring and integrating over the time 
interval T which is given by [l26l eq. (2)] as 



( 2 ) 
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The output of the integrator is a measure of the energy of the reeeived waveform whieh 
eonstitutes a test statistie that determines whether the reeeived energy measure eorresponds only 
to the energy of noise (Hq) or to the energy of both the unknown deterministie signal and noise 
(Hi). By denoting the observation time bandwidth produet as u = TW, the test statistie follows 
the eentral ehi-square distribution with 2u degrees of freedom under the Hq hypothesis and the 
non eentral chi-square distribution with 2u degrees of freedom under the Hi hypothesis H]. To 
this effect, the corresponding probability density function (PDF) in the presence of AWGN is 
expressed according to eq. (3)], namely, 


where 
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is the corresponding instantaneous SNR and Eg denotes the signal energy. Furthermore, 


( 4 ) 


r (a) = / 


( 5 ) 


stands for the Euler’s gamma function and 


4(a:)4l/' cos(n0)e"“^(®)d0, (6) 

^ Jo 

denotes the modified Bessel function of the first kind ll49l . 

As already mentioned, an energy detector is largely characterized by a predefined energy 
threshold, A. This threshold is particularly critical in the decision process and is promptly asso¬ 
ciated to three measures that overall evaluate the performance of the detector: i) the probability 
of false alarm, Pf = Pr(r/ > A | Hq); ii) the probability of detection, Pd = Pr(r/ > A | Hi) and 
iii) the probability of missed detection. Pm = 1 — Pd- The first two measures are deduced by 
integrating ([3]) over the interval between the energy threshold to infinity, {A, cxo}, yielding 

p r (u, I) 

' r(«) ■ 


( 7 ) 
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and 


Pd = Qu{y/^, \/A), 


( 8 ) 


where, 


f 

r {a,x) = / 


oo 



(9) 


denotes the upper ineomplete gamma funetion ll^ and 



( 10 ) 


is the generalized Mareum Q—function Il50l . 

B. The Weibull Fading Distribution 

As already mentioned, the Weibull fading model has been shown to represent effectively 
the small-scale variations of a signal in non-line-of-sight (NLOS) communications in multipath 
environments. Physically, this fading model accounts for propagation in non-homogeneous en¬ 
vironments and the resulting envelope is obtained as a nonlinear function of the modulus of the 
sum of the involved multipath components [|4^ . The envelope PDF of Weibull distribution is 
given by ll4^ eq. (2.27)], namely. 



( 11 ) 


where a denotes the fading parameter that is selected to yield the best fit to measurement results 
and D = = £'[/?“] is the mean power of the desired signal H^ . The corresponding PDF of 

the instantaneous SNR per symbol 7 is given by [l43l eq. (2.29)], namely. 



( 12 ) 


where 7 represents the average SNR per symbol. To this effect, the corresponding envelope and 
power CDFs are expressed as. 



( 13 ) 
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and 


P^(7) = 1 - 


respectively Il43]l . The statistical moment is expressed as, 


E[R^] = r 


1 + 
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aJ 


whereas the corresponding MGF is given by [|^ eq. (2.35)], 
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where G;|;(.) denotes the Meijer G—function If54ll . 
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III. Average Detection Probability over Weibull Fading Channels 

It is recalled that (|7|) and ([ 8 ]) account for the case of AWGN channels. For communication 
scenarios over fading channels the average probability of detection is obtained by averaging ([ 8 ]) 
over the corresponding SNR fading statistics, namely. 


(17) 

Based on this, the average probability of detection in the case of Weibull fading channels can 
be obtained by averaging ([ 8 ]) over the statistics of (fT^ . To this effect, by substituting (fT^ in 
(fTTl) it immediately follows that. 
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Evidently, an analytic expression for the Pd in Weibull fading conditions is subject to analytical 
solution of the integral in (fT^ which is not available in tabulated form in the open technical 
literature. To this end, we firstly let x = 1/27 and thus: 7 = x^/2 and dx/d'j = By 

substituting in (fTSl) and after some basic algebraic manipulations one obtains. 
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By integrating Xi by part it follows that, 
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By recalling that Qn(oo, 6 ) = 1 and based on [|5T1 eq. ( 2 )] and [[STl eq. (10)], it immediately 
follows that. 


<340. ^/A) = e -4 x; X (22) 
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and 
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respectively. Thus, by substituting (l22l) and (1231) in (|2T]) yields. 
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By expanding the exponential term with the arbitrary power, the X 2 term can be equivalently 
written as follows. 
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Importantly, the X 3 integral can be expressed in closed-form according to lEH eq. (6.643.2)]. To 
this effect, by making the necessary change of variables and substituting in (|25]) and subsequently 



























in (l24l) and (fT9l) . the following analytic expression is deduced for the average probability of 
detection, 
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is the Kummer’s confluent hypergeometric function and 


r(a + n) 

r(a) 


(28) 


denotes the Pochhammer symbol I1491I . It is noted that the series converges rapidly and only few 
terms are required for achieving acceptably low levels of truncation error. This is particularly 
important since spectrum sensing applications typically require accuracy of two decimal digits. 

To the best of the Authors’ knowledge, equation (|2^ has not been reported in the open 
technical literature. 


IV. Numerical Results 

Having derived a novel analytic expression for the average probability of detection, this section 
is devoted to the analysis of the behaviour of energy detection over Weibull fading channels. 
The corresponding performance is evaluated for different scenarios of interest through both Pd 
versus 7 curves and complementary receiver operating characteristics (ROC) curves {Pm versus 
Pf). To this end. Fig. 1 demonstrates Pd vs 7 curves for different values of fading severity a 
with Pf = 0.1 and u = 5. One can clearly notice that the behaviour of the energy detector is 
very sensitive as even at slight variations of a, the value of Pd changes rapidly. As expected, 
the performance of the detector is directly proportional to the value of the severity index a and 
the average SNR 7. For small values of a and 7 the performance of the detector for Pf = 0.1 
is rather low and thus, increasing the time-bandwidth product u is essential. Likewise, Fig. 2 
depicts the behaviour of Pd vs 7 for higher values of a. The performance of the detector increases 
significantly, particularly in the post-Rayleigh area i.e. a > 2. Furthermore, the overall flexibility 
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Fig. 1. Pd vs 7 for i.i.d Weibull fading with P/ = 0.1, u = b and different values for fading severity a. 
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Fig. 2. Pd vs 7 for i.i.d Weibull fading with P/ = 0.1, u = b and different values for fading severity a. 
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Fig. 3. Complementary ROC curves for u = 5, 7 = —5dB and different values of a. 


of the Weibull model is evident sinee it can account for a wide range of severity conditions that 
have a significant effect on the corresponding performance. 

Figs. 3—5 illustrate complementary receiver-operating-characteristics (ROC) curves vs Pf 
for M = 5 at: i) low SNR (7 = —5dB); ii) moderate SNR (7 = lOdB); in) high SNR (7 = 25dB) 
and different severe fading conditions. It is clearly shown that the effect of severe fading on 
the system performance is stronger as the average SNR increases. For example, for Pf = 0.2 
and a = 1.0, the corresponding value for the probability of missed detection is Pm = 0.78, 
Pm = 0.41 and Pm = 0.1, respectively. However, it generally appears that the detector can 
perform adequately even at significantly more severe fading conditions that Rayleigh fading. 
This also verifies the usefulness of Weibull distribution as an adequate model for accounting for 
the multipath fading effect. 


V. Conclusion 


This work analyzed the performance of energy detection over Weibull fading channels. A 
novel analytic expression was derived for the average probability of detection and was shown 
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Fig. 4. Complementary ROC curves for u = 5, 7 = 10 dB and different values of a. 



Fig. 5 . Complementary ROC curves for u = 5 , 7 = 25 dB and different values of a. 
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that the overall performance of the detector is largely affected by the value of the involved 
parameters. This is evident by the fact that even slight variations of the fading parameter results 
to substantial variation of the corresponding probability of detection. As a result, the offered 
results are useful in quantifying the effect of fading in energy detection spectrum sensing which 
can contribute towards improved and energy efficient cognitive radio wireless communication 
deployments. 
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